Phase Transitions in the One-Dimensional Pair-Hopping Model: 
a Renormalization Group Study 
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The phase diagram of a one-dimensional tight-binding model with a pair-hopping term (amplitude 
V) has been the subject of some controvery. Using two-loop renormalization group equations and 
the density matrix renormalization group with lengths L < 60, we argue that no spin-gap transition 
occurs at half-filling for positive V, contrary to recent claims. However, we point out that away 
from half-filling, a phase-separation transition occurs at finite V. This transition and the spin-gap 
transition occuring at half-filling and negative V are analyzed numerically. 



PACS Numbers: 75.20.Hr 



I. INTRODUCTION 



In recent years interest in correlated electron systems 
has increased particularly in an attempt to understand 
high-temperature superconductors. It is important to 
study the effect of all possible nearest-neighbour interac- 
tions; however, the hopping of on-site spin-singlet pairs 
has not been well studied to date. 

The Hamiltonian for the pair-hopping model is Q 



H 



h.c] 



<«■> 



(i) 



where a- (ci CT ) creates (destroys) an electron of spin a = [ 
, i at lattice site i, so that rii a — c\ a ci a is the number of 
electrons of spin a at site i, and (ij) denotes nearest- 
neighbour pairs. Thus t and V are the single-electron- 
and pair-hopping amplitudes respectively, so this models 
a competition between the two hopping terms. As t — ► — t 
is a symmetry of H, in this paper we take t > 0. 

At large \V\/t, all sites are doubly occupied or empty 
(assuming an even number of electrons) and the model 
becomes equivalent to spinless fermions jl] . In particular 
there is a large gap, of 0(V), to any excited state with 
non-zero spin. This is true for either sign of V, but it 
is important to note that V — > — V is not a symmetry 
of the model, unlike the Hubbard model. Finite-size nu- 
merical work has been performed on this model in 
one dimension for positive V, suggesting a phase transi- 
tion at which the spin gap (or single-particle excitation 
gap) opens, at V/t w 1.4. Two different analytical renor- 
malization group (RG) analyses have been applied to the 
model. One (2| suggested the existence of a spin gap A s 
for all V > 0, with 



oc e 



-7Tt/V 



(2) 



as V/t — > 0, and no transition for any positive V; the 
other H suggested that there is a transition at V/t w 1.4, 
consistent with the numerical work. 
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One purpose of this paper is to re-examine this ques- 
tion. Previous numerical work has used chains of length 
L < 12. We present data for much longer chains, 
L < 60, using White's density matrix renormalization 
group (DMRG) technique thereby countering the 
dominance of finite-size effects. We also discuss the sub- 
tleties involved in trying to extract information about the 
phase diagram from a low-order analytical RG calcula- 
tion. Our conclusion is the same as that of Ref . 0j : no 
phase transition for any V > 0. 

It was argued in Ref. 0] that there is a phase transi- 
tion, corresponding to the appearance of a spin gap, for 
some finite negative V. (The case V < was not studied 
in Refs. jl| or jj|.) We analyze this numerically, finding 
a transition at V c rj —1.5t. 

There has been considerable interest of late in phase 
separation in the Hubbard and t-J models in one and 
higher dimensions. The pair-hopping model provides a 
simple example of a model where it is easy to see that 
a phase-separation transition must occur at some finite 
critical coupling, with a non-zero total magnetization. 
Consider the model at t = 0, with a total magnetization 
M/2, corresponding to an excess of M spin-up electrons. 
These electrons necessarily reside on singly occupied sites 
and are therefore completely immobile when t = 0. Thus 
the model is equivalent to spinless fermions with vanish- 
ing hopping terms to M sites. Equivalently we have an 
XY spin chain with vanishing exchange coupling to M 
sites, corresponding to M nonmagnetic impurities. These 
simply have the effect of breaking the chain up into chain- 
lets. It is fairly clear, and can be demonstrated explicitly 
by a trivial calculation in the free spinless fermion model, 
that the energy is lowest when all M impurities are next 
to each other, leaving an XY chain of L — M sites. This 
corresponds to phase separation at \V\/t ^> 1 for any 
non-zero magnetization. On the other hand, the renor- 
malization group analysis of the model at weak coupling, 
\V\/t <C 1, indicates behaviour similar to that of the Hub- 
bard model, with no phase separation. This suggests that 
phase-separation transitions should occur at finite values 



1 



of V/t (one for positive V and one for negative V). We 
find evidence for such a transition at V = V c i ~ 3.5f, but 
we have not examined the one at V = V C 2 < 0. 

In Sec. [il] we review and critically analyze the previ- 
ous analytical RG calculations. In Sec. Ill we analyze 
quantitatively the large \V\/t limit. In Sec 



IV B 



we 



present our numerical work at half-filling, M = and 



V > 0, indicating no phase transition. In Sec. IV C we 
present numerical evidence for the spin-gap transition at 
half-filling, M = and V < 0. In Sec. pVD| we present 
numerical evidence for the phase-separation transition for 
M ^ at V = V cl w 3.5t 



II. ANALYTICAL RG STUDIES OF THE PHASE 
DIAGRAM 

The RG analyses of Refs. Q and || came to quite 
different conclusions. Here we would like to explain the 
reasons for this and give arguments in favour of the for- 
mer approach. 

We use essentially the notation of Ref. S, which is 
taken from the review article of Solyom la]. Taking the 
continuum limit of the pair-hopping model we obtain a 
general Hamiltonian: 



h = J dx[H + n- mt ], 



(3) 



where Ho and Hint are the dimensionless kinetic energy 
density and interaction Hamiltonian density. We keep 
only wave vectors close to the two Fermi points, k « ±fcp, 
which we label, in position space, ipi, and ^r,. The kinetic 
energy density is given by 

Ha = v F [tl)"h—ijj La - V ) R t i^V , Ra]- (4) 

Here the spin index a is implicitly summed over; v F is 
the Fermi velocity. Some of the various interaction terms 
can be conveniently written in terms of charge and spin 
currents (or densities) 



J[ ee and Jl ee Vl T ^l,3, 
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(•5) 



and similarly for and Jj|. In this notation, we have 



[e Q/ #£ T < T e 7 VR 7 ^ + (L <- R)] 



We have chosen to write the Hamiltonian in a manifestly 
SU(2) invariant way. The last term can also be written 
as 



« - 3^ ' JL = KJl ~ 4«* = 4Jl T Ju, (7) 
where 

Jlq = V'l^Lo (repeated index not summed). (8) 

To the first non- vanishing order in V, the bare couplings 
have the values 

v F = 2t; g p = -g s = g 3 = g 4 = 2V/ttv f . (9) 

To cubic order, the RG equations are given by 

-^=S S 2 + K<?s + <?4)<? S 2 

~~df = 9p93 + 3(5p + ffl - 2g p gi)g z 
ji — 4W93 y s )- 



(10) 



Here I = — log A, where A is an ultraviolet cut-off. As we 
lower the cut-off to study the long-distance behaviour, I 
increases. 

Part of the discrepancy between the conclusions of 
Ref. [0] and Ref. || arises from the treatment of the g 4 
coupling. If we bosonize the theory, then 



Z7T 



4J£ Z J£ Z 



2tt 



(11) 



where PjS are charge and spin bosons. Hence 174 simply 
shifts the velocities of charge and spin excitations to 



v p = «f(1 +34/2) 
v s = v F (l - 54/2). 



(12) 



A common approach to Luttinger liquids is to sim- 
ply set v p and v s to their renormalized values and drop 
174 from the RG equations. This approach was used in 
Ref. p|. The RG equation for g s then decouples from 
the g p and g 3 ones. This arises from the fact that, upon 
bosonizing, the corresponding operators involve only the 
spin boson and only the charge boson respectively. We 
then see that g s = is not a stable fixed point: if g s < 0, 
as is the case for V > 0, g s will flow away to strong cou- 
pling. This is usually taken to indicate that the system 
is in a phase with a gap for spin excitations. 

On the other hand, a quite different conclusion can be 
reached if g 4 is kept in the RG equations. Then, accord- 
ing to Eq. ( |To| ) , g s — becomes a stable fixed point from 
the negative side provided that g 4 < —2. The nature 
of this putative phase can be understood by also rewrit- 
ing the free electron kinetic energy in terms of spin and 
charge currents. Setting all coupling constants to zero 
except 174, the full Hamiltonian density can be written 
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n 



TTV-p 



(i-f)« + (i + f) 



J T s v T s 



(L 



(13) 



We see that for g 4 < —2, the spin part of the Hamilto- 
nian becomes unstable. That is, J£ 2 (a;) and J§^(x) tend 
to become large, necessitating the keeping of higher or- 
der terms in the Hamiltonian. On the other hand, the 
condition of zero total magnetization requires 



dx[Jl 



Jr] = o. 



(14) 



A possible interpretation of this phase (which occurs in 
other known cases) is a ferromagnetic phase. The con- 
dition of zero total magnetization forces a domain struc- 
ture, i.e. phase separation, to occur. One side of the 
system has positive polarization and the other half neg- 
ative. 

In Rcf. S, the cubic RG equations were integrated, 
including g 4 , using the initial values of Eq. (0) [plus the 
0(V 2 ) corrections which are not important at small V]. 
The result was that fui < Vft^^b, a fiAcd puinL was 

reached with g 4 « —2.5 and g s = 0. 

Whether or not 34 is included, for V/t < 1, g% renor- 
malizes to zero and g p to some small positive value which 
depends on V/t, corresponding to zero gap for charge ex- 
citations. 

In Ref. Q this phase was identified as having a spin gap 
since g s does not flow to zero. In Ref. || this phase, with 
g s = and g 4 < —2, was assumed to have no gap for 
single-particle excitations. Since these excitations have 
spin i and charge 1 this would imply, from the usual Lut- 
tinger liquid viewpoint that there is neither a charge gap 
nor a spin gap. We do not find this calcuation convincing. 
It is not possible to argue rigorously that 54 renormal- 
izes to a value less than —2 using only the cubic order 
RG equations. If this actually happened, as claimed in 
Ref. [|3| , this would presumably imply a transition into a 
ferromagnetic phase (or possibly some other more exotic 
phase characterized by the harmonic spin Hamiltonian of 
Eq. (|l^) becoming unstable) for arbitrarily small V. No 
direct numerical evidence for ferromagnetism (or other 
exotic behaviour) at small V/t has been presented. Al- 
though earlier numerical work in Refs. p]] and H saw 
indications of vanishing spin gap in this region of param- 
eters, the numerical results presented here in Sec. IV B 
based on much longer chains (L < 60 instead of L < 12) 
find a non-zero spin gap. 

In Ref. H a different phase is reached for V/t > 1 with 
a non-zero g s at the fixed point, corresponding to a spin 
gap as in Ref. O. However Refs. ^ and || now disagree 
about the behaviour of the charge couplings, g p and g%. 
Note that the second and third RG equations in Eq. jlO| ) 
imply that 53 = ±g p are separatrices (for g 4 = 0). For 
g p > 0, if |.93 1 < g p , g 3 flows to zero (see Fig. [l]), cor- 
responding to a harmonic gapless effective Hamiltonian 




(b) 
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FIG. 1. Third-order RG flow diagrams, ignoring g 4 . 

(a) Flow in the charge sector, for small gs and g p . 

(b) Flow in the spin sector. 



for charge. Outside this region both 33 and g p flow off to 
values of 0(1). This is normally interpreted as a phase 
with a charge gap. It is a remarkable feature of the pair- 
hopping model that, to 0(V), g p = .93: the system lies 
on a separatrix. It is necessary to calculate the bare 
couplings to 0(V 2 ) to deduce whether or not 173 flows 
to zero. Both papers agree that these 0(V 2 ) terms place 
the bare couplings in the basin of attraction of the gs = 
critical line, for small V. In Ref. it was assumed (on 
the grounds of simplicity) that the system remained in 
this basin of attraction for all V > 0. On the other 
hand, in Ref. [|), the expression for the bare couplings 
to 0(V 2 ) was used for arbitrarily large V to deduce that 
the bare couplings moved outside this basin of attraction 
at a critical V ~ t (the same critical point at which g s 
and g 4 change). The cubic RG equations predict a fixed 
point at g p = 93 = —2, which the authors of Ref. Q 
assume corresponds to vanishing charge gap. 

This argument concerns us because it is not possible 
to tell from these low order calculations of the bare cou- 
plings and the RG equations whether or not the bare cou- 
plings ever leave the domain of attraction of the 93 = 
critical line. Furthermore, if they did, this phase would 
normally be identified as having a charge gap, which we 
know does not occur for small or large V. (The existence 
of an apparent, finite coupling fixed point of the cubic 
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RG equations at couplings of 0(1) does not necessarily 
signal the existence of a different critical point. It could 
disappear upon keeping higher order terms.) 

In Sec. Ill we give analytic arguments implying that, 
for large \V\/t, there is a spin gap but no charge gap. 

By ignoring g^ (i.e. absorbing it into velocity renor- 
malizations) and making a plausible assumption about 
the behaviour of bare coupling constants at large V, we 
obtain simple behaviour requiring no phase transition for 
any V > 0. There is always a spin gap and no charge 
gap. 

On the other hand J^] , by including the renormalization 
of 174 and using weak coupling results at strong coupling 
one obtains two different phases: a bizarre small V phase 
with an unstable harmonic spin Hamiltonian and a large 
V phase which would likely correspond to a charge gap, 
in contradiction with the expected large V result. 

The authors of Ref. H applied the same RG analysis 
to the positive U Hubbard model. Their analysis gave a 
small U phase with 94 > 2, corresponding to a negative 
harmonic Hamiltonian in the charge sector and a large 
U phase with a non-zero g s which would normally corre- 
spond to a spin gap. As they pointed out themselves, this 
is in contradiction with the expected behaviour which is 
a charge gap and no spin gap for all positive U (at half- 
filling). 



III. LARGE \V\ LIMIT: SPIN GAP AND PHASE 
SEPARATION 



The pair destruction operators a, = Cj^Cif have com- 
mutation relations [aj,at] = (1 — riijSij, where n, = 
^2 a n io-- So the a's commute if rii = 1 since they have 
no effect on singly occupied sites. In the other two cases 
n-i = 0,2 we have n% = 2a\a,i so that the a^'s are spinless 
fermion operators, obeying {ai,aj} = Sij. As shown in 
Ref. (I), setting t = results in a ground state involv- 
ing only empty and doubly occupied sites, so the on-site 
pairs are effectively spinless fermions. 

The DMRG method that we use requires free bound- 
ary conditions. So we analytically examine the large \V\ 
limit for an open chain, noting that nothing essential will 
change in going to periodic or infinite length chains. The 
ground state energy for an open chain of even length L 
is easily computed to be 



L/2 

E Q = -2\V\ cos T~~~T = l^ 1 ( 1 ~ csc T 

n=l 



2L + 2 



(15) 



Adding a single electron to this half-filling ground state 
produces an immobile site since t — 0, effectively break- 
ing the chain. The energy will depend on the location 
of the break, and is easily shown to be minimized if the 



break is at the end of the chain, in which case the energy 
is that of h pairs hopping on an open chain of length 
L — 1, namely 



L/2 



E 1 = -2|y|^cos^ = H (l 



cot • 



2L 



So the single-particle gap for the open chain is 



Asp = \V\ 



W\ 



CSC • 



2L^ 

2 TT 



TT TT 

cot 

2 2L 



TT 



12L 2 



(16) 

(17) 
(18) 



So for t = 0, we have a model equivalent to free spinless 
fermions, corresponding to a spin gap proportional to 
|V| but no charge gap. To see whether this situation 
persists for finite |V|/t, we can do perturbation theory 
in the lattice model in t/V. This is very similar to the 
well-known results on the large-[/ Hubbard model. In 
this case we project out singly occupied sites. A single 
application of t takes us into the high-energy subspace 
with 2 singly occupied sites. In second order perturbation 
theory we generate an effective interaction of 0(t 2 /V) in 
the spinless fermion model. This simply corresponds to 
a nearest-neighbour interaction of the spinless fermions. 
This interaction is known to be exactly marginal, leading 
to a critical line with vanishing gap. 

Thus there is a spin gap for large |V|/t. As there is 
no spin (or charge) gap for V = 0, there must be some 
transition. On the basis of a reliable interpretation of 
the analytical RG equations and careful consideration of 
and comparison with numerical RG results, we conclude 
that the positive V transition occurs at V = instead of 
at some finite \V\. We show numerically that for small 
positive V, the behaviour of the single-particle gap is of 
the form predicted by the RG flows (upon dropping g±) 
in the numerically accessible region of phase space. We 
also find a spin-gap transition at V — V c w — 1.5t. 

The above analysis also shows that in the case t = 0, 
a single unpaired electron sits at a chain end; it is clear 
that additional electrons of the same S z will clump at 
the chain ends as well. That is, at finite magnetization, 
the chain phase separates: one part of the chain assumes 
the net magnetization. It is important to note that this 
is not a peculiarity of the open chain; in the periodic 
case as well, at t — 0, added polarized electrons cut the 
chain and the chain-breaking energy is clearly minimized 
by clumping them together. Since going from t = to 
some large but finite \V\/t introduces only a marginal 
operator, it is clear that this phase separation will per- 
sist to some critical values of V, which probably have 
different absolute values because the Hamiltonian is not 
symmetric under V — > — V. 
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IV. RESULTS OF NUMERICAL RG 



A. DMRG Details 

We use the "infinite system DMRG method" @] , treat- 
ing open chains of even length up to 60, and maintain- 
ing 64 (sometimes 128) states in each block. The ground 
state has total spin and is at half-filling; we add a single 
electron (pair of electrons of opposite spin) to compute 
the spin (charge) gap for each length. These results are 
extrapolated to infinite length taking into account trun- 
cation error uncertainties. The figures summarize the 
results of our DMRG calculations, as explained in this 
section. 



B. Spin Gap for V > 

We find that the spin gap does not vanish for any V > 
0, as shown in Fig. |2|. In comparing its dependence on V 
with that predicted from the analytical RG of Ref. @, 
namely 



te 



-■Kt/V 



(19) 



(after correcting the typographical error) which is valid 
for small V/t, the numerical work is not dependable for 
V/t < 1 because there the expected correlation length 
£ w v-p/2A becomes of order the system size L. The 
finite-size gap alone is Aps ~ itv-f/L so that one cannot 
expect to measure A/t lower than Aps/i « 2ir/L « 0.1 
for L = 60. 

The RG flow equations to two-loop order, after drop- 
ping (74 as explained in Sec. [n], give for the spin coupling 



9s0 



— log - 

2 l + 29 s - 



h - = \og — 

^0 



(20) 



where Lq is an initial length scale (A = L 1 is the ultra- 
violet cut-off) and the initial spin coupling is || 



V 



V\\ ( 1 

— log tan 

ntj & \ tL 



(21) 



We take the spin gap to be the energy (inverse length) 
scale at which g s enters the regime of strong coupling, 
specifically where g s = a = 0{— 1), resulting in 



A S (V) = A exp( ff: 



sQ 1 




(22) 



where a = 0{— 1) is used as the criterion for \g s \ becom- 
ing large. Fig. || shows that the numerically computed 
spin gap is indeed of the form predicted by the RG flow 
equations with 34 dropped. 




FIG. 2. Summary of numerical results: the open squares 
and dashed line are the charge gap, and the filled squares 
and solid line are the spin gap. A clear phase transition is 
evident near V — V c ~ — 1.5t, but for positive V, the spin gap 
opens up from V = 0. The error bars indicate uncertainty in 
extrapolating L^ 1 — > 0; the lines are to guide the eye. 
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-2 -1.5 -1 -0.5 

-t/V 

FIG. 3. Fitting the DMRG data for only two points, 
namely V/t = 1, 2, to the form given by Eq. (p^), the dashed 
lines are the upper and lower limits of the resulting fitted 
curves taking the numerical error bars into account. The 
lower limit extrapolates well over the range t < V < 4t, which 
is the expected region of validity. (The fit is not expected 
to be valid for V < t because of the finite-size gap, while 
V > Vd ~ 3.5t is the phase-separated region.) 
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While the above comparison of the RG flows are to 
numerical results on open (not periodic) chains, we be- 
lieve these results (and in particular the non-vanishing 
of the gap for V > 0) constitute a reliable estimate of 
the situation in the thermodynamic limit. However the 
situation is somewhat different in the phase-separated 
region (V > V c i « 3.5i and V < V c2 < 0) than in the 
non-phase-separated region at smaller \V\. In the non- 
phase-separated region, the excitation which we study 
is concentrated in the bulk of the chain as discussed in 
Sec. 1VD and Figs. H-El Thus we expect that its exci- 



0.25 



V/t 



tation energy is not affected significantly by the bound- 
ary conditions for sufficiently long chains. However, in 
the phase-separated region, the excitation lives near the 
ends of the chain and its energy may well be strongly 
affected by the boundary conditions. In this case, the 
energy which we measure is still a lower bound on the 
bulk gap. This follows because the state which we study 
is the lowest energy one with these quantum numbers. If 
the bulk gap were lower, we would expect a lower energy 
state to exist, localized far from the chain ends. Thus our 
results give strong evidence for a spin gap for all V > 
but only give a reliable estimate of the size of the gap for 
V C 2 < V < V c i ~ 3.5i, except for magnitudes less than 
the finite-size gap as discussed above. 



C. Phase Transition at V — V c 



-1.5* 



As discussed in Rcf. 0], for small V < the pair- 
hopping model is identical to the positive U Hubbard 
model. Thus we expect a charge gap but no spin gap in 
this region. It was also argued that there should be a 
phase transition at finite V < because at V/t — > — oo 
there is no charge gap but a spin gap. In Fig. || we 
present DMRG results confirming this prediction, with 
the transition occurring at V = V c « — 1.5t. Our numer- 
ical results are consistent with the spin gap appearing 
at the same critical coupling at which the charge gap 
disappears; however, the presence of two distinct critical 
couplings cannot be ruled out. It is unclear to us whether 
this critical point (or points) simply corresponds to the 
renormalized couplings g s and 93 passing through zero, 
or to some more exotic critical point. 



D. Phase Separation at V = V c i ~ 3.5t 

To demonstrate the phase-separation transition, we ex- 
amine the behaviour of wave functions obtained using the 
DMRG at L = 60 in the sector of one electron added rel- 
ative to half-filling. Specifically, we plot in Fig. ^ the ex- 
pectation value of S z (i) for sites i = 1, . . . , 30 (the chain 
is symmetric about the central link) for different values 
of V/t. For large V, the excess spin is localized at the 
chain ends, and as V is reduced, the spin extends further 



V 




20 

Site Index i 

FIG. 4. Expectation values (S z (i)} for for different values 
of V/t, for one electron added relative to half-filling. The 
unpaired electron delocalizes into the chain near V/t — 3.5. 
(The L — 60 chain is symmetric about its central link.) 



into the bulk. As V/t drops from 4 to 3, looking at the 
wave function near the centre of the chain shows that 
near these values of V/ 1 the spin becomes unbound from 
the chain end and is rapidly and fully delocalized into 
the bulk of the chain, leading us to consider V c \ « 3.5t 
as a phase-separation critical point. This conclusion is 
further verified by examining the spin on the chain end 
as a function of V/t, as well as the total spin in the centre 
half of the chain, as shown in Figs. || and 

Due to the fact that we have employed the infinite 
system method, instead of the finite system method [Q, 
these wavefunctions are not expected to be precise par- 
ticularly near the phase-separation transition and at the 
chain ends. However, we expect that the results are ac- 
curate to within a few per cent at worst, certainly not 
affecting the qualitative behaviour of our figures which 
clearly demonstrate the phase-separation transition. 

While a priori this phase transition could occur at a 
different value of V than the bulk phase separation, the 
simplest scenario would have both transitions occurring 
at the same point: essentially the bulk transition drives 
the boundary transition. The numerical evidence on one 
and two added electrons seems to indicate that for low 
net magnetization, V c \ is constant. 

This phase-separation transition will occur for finite 
V c i in the periodic and infinite chain as well (though 
not necessarily at the same value of V c \ as for the open 
chain) : added unpaired electrons will still break the chain 
into chainlets and the energy will be minimized if they 
clump together. However, it will be more difficult to 



G 



detect in a periodic chain since the ground state is usually 
translationally invariant. 




2 4 6 8 



V/t 

FIG. 5. Spin at a chain end (open squares) and net spin in 
the centre half of the chain (filled squares) as a function of 
V/t for a single added electron. 



V. CONCLUSIONS 

We conclude that there is a finite spin gap for all pos- 
itive V in the half-filled pair-hopping model in one di- 
mension, and that to accurately describe its behaviour 
as a function of V, one must neglect the coupling 174 in 
the renormalization group flows. 

We conclude that there are phase-separation transi- 
tions in the pair-hopping model, one at positive V and 
one at negative V. In one dimension at low doping from 
half-filling, for V > V c i ~ 3.5i polarized electrons clump 
together. 

We conclude that there is a new critical point at V = 
V c « — 1.5t at which, proceeding from weak coupling, a 
spin gap opens and the charge gap closes at half-filling. 
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2 4 6 8 

V/t 

FIG. 6. Spin at a chain end (open squares) and net spin in 
the centre half of the chain (filled squares) as a function of 
V/t for two added electrons. 
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